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INTEGRABLE Z„-CHIRAL POTTS MODEL: 
THE MISSING RAPIDITY-MOMENTUM RELATION 



G. VON GEHLEN " 
Physikalisches Institut der Universitdt Bonn, 
Nussallee 12, D-53115 Bonn, Germany 

The McCoy-Roan integral representation for gaps of the integrable Z„- 
symmetric Chiral Potts quantum chain is used to calculate the boundary of 
the incommensurable phase for various n. In the limit n ^ oo an analytic 
formula for this phase boundary is obtained. The McCoy-Roan formula 
gives the gaps in terms of a rapidity. For the lowest gap we conjecture the 
relation of this rapidity to the physical momentum in the high-temperature 
limit using symmetry properties and comparing the McCoy-Roan formula 
to high-temperature expansions and finite-size data. 



1 Introduction: Z„-symmetrical Chiral Potts Quantum Chains 

The Chiral Potts model has been introduced in 1981 by OstlundB in order to describe 
commensurate-incommensurate (C-IC) phase transitions observed in surface monolayer 
adsorbates, e.g. krypton on a graphite surface. The location of a possible Lifshitz point 
in the phase diagram was a maior issue. The quantum chain version has been used 
first by Centen et azi. After the discovery of the novel integrability properties of 
the Chiral Potts modelsHi, much interest has focussed on their mathematical aspects. 
Nevertheless, the spectrum of the Chiral Potts quantum chains exhibits several unusual 
interesting physical features caused by the inherent parity violation. 
The Z„-symmetrical Chiral Potts quantum chains are defined by 

N n-l 

Hi-) = _ ^ ^ {a, + A a, rf r,^t } . (1) 

i=l k=l 

The operators cjj and Fj act in a vector space C" located at site i and satisfy 

a^,=rja.J^^^; a: = r^ = l; 4 = a^-\ r+ = (2) 

where oj = e^'^*/". The parameter A > plays the role of an inverse temperature (we 
shall consider the ferromagnetic case only), and are complex parameters. If we 
choose to represent the by diagonal matrices {ai)i^m = ^''~^^i,m then the Fj become 
cyclic lowering operators (rj)/ ^ = (J/ ^+i- is the number of sites. We use periodic 
boundary conditions Tn+i = Ti. 

A three- parameter version of (||), expressing the coefficients in terms of just two 
chiral angles < 99, < vr: 

afc = e*'^(^~^Vsin(7rA;/n); Ofc = e^'^^^'^V sin (7rA:/n) (3) 

describes various cases of particular interest: 
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For general n and ip = (p = we recover the Fateev-Zamolodchikov parafermionic 
model, for n = 2 we get the standard Ising quantum chain. 



li if = (p = 7r/2 we get a representation of the Onsager-Dolan-Grady-algebra 



Lt (f = (p = IT/ / we get a representation 
"Superintegrable Chiral Potts-Model" ii. 



• An- Yang et al. i have shown that (||),@ can be derived from a two-dimensional 
lattice model with (generally complex) Boltzmann weights satisfying Yang-Baxter 
equations (so that H^'^'^ is integrable) if the three parameters are related by 

cos (/3 = A cos (p. (4) 

• Eq.(@) forces to become non-hermitiar^ for A < 1 if we choose | cos ip\ > 
A. For (|^ with cos if > X the hamiltonian is parity invariant up to complex 
conjugation, whereas for (|^ with cos (f < —X there is a more complicated shifted 
parity symmetry see eq.(15) below. 

H^^') is Z„-symmetrical and translational invariant. So, its spectrum decomposes into 
sectors labeled by the Z„-charge Q = 0, . . . , n — 1 and P, the momentum eigenvalue 
which takes the values < P = 2Trk/N < 27r with A; = 0, 1, . . . , - 1 or -[A/2] <k< 
[N/2]. 

We denote the eigenvalues of i/^") at fixed A, ip (in this talk (p will always be 
determined by (§)) by Eq^riP)- The index r = 0, 1, . . . labels the levels in each Q, P- 
sector, starting with r = as the ground-state of the respective sector. In some part, 
but not in all regions of the A, 99-plane, the ground state of H is in the Q = P = 0-sector. 
By convention, we shall define gaps always with respect to the lowest Q = P = 0-level 
at the same A, (p and ip: 

Ai?Q,,(P) =: EQ,r{P) - Eo,o{P = 0). (5) 

Since in general parity is not conserved, a relation E{P, A, ip) = E{—P, A, ip) is expected 
to be valid only for special values of A and ip. Consequently, in general the minimum 
of E{P, A, ip) is not at P = as is seen e.g.^iiL_Eig. 2. Quite generally, the spectrum of 
seems to be a quasiparticle spectrum with one quasiparticle in each Q ^ 0- 

charge sector and one or more exceptional "Cooper-excitations" @. In the IC phase 
the ground state looses translational invariance, because there the Q = 1-quasiparticle 
dispersion curve dips below the lowest Q = P = 0-level for a certain range of P 7^ 
and AEQ=i^r=oiP) becomes negative. 



2 The McCoy-Roan gap formula 

While for the Za-superintegrable case the complete spectrum of H^"^^ has been calcu- 
lated analyticallyQ'cl, only few details are known about the spectrum in the more general 
integrable case (^). In the rest of this talk, we shall consider this integrable case. The 
difficulty of making progress here comes because in this case the rapidity parameters 
appearing in the Yang-Baxter-equations are restricted to higher elliptic curves. Baxter 
t3 has devised an approach which avoids uniformization of the spectral variables. Solv- 
ing a functional equatioilp due to Baxter, Bazhanov and PerktJ, McCoy and Roan for 
Z3 , and then McCoy for the general Z„-caselll have derived an analytic formula for the 



low-temperature gap in the thermodynamic hmit. By duahty, this is equivalent to the 
following expression for the lowest AQ = 1 gap in the high-temperature regime A < 1: 



A-^-l){(l + A)2-(l-A)2t«} 



V /-(i^x)'^; 2(tt;-cosf) V(i"-l){(l + A)^-(l-A)^t«} 



- 2t7;cos^ - 1 



(6) 



Here P in front of the integral indicates the principal value. The sign it in front of the 
second term has to be chosen " -|- " for < (/> < ^ and " — " for ^ < < vr. The chiral 
angle (j) enters @ only in the combination 

/(A,(/)) = l-2Acos((/>-99) + A^ (7) 

through tp = /(A, </>)/(! — A)^. is a rapidity parameter which will also be used in the 
form V = {{1 + A)/(l - A))^/"^;. The relation of 11 or t; to the physical morrientum P 
will be our main concern. Only in the superintegrable case it is known to read 
exp(-iP) = {l + uj'^v)/{l +ujv). 

We find it useful to recast into the following form: 



AEi^oiv, A) = 2 J/(A, 0) 1 + 2 Qm- 



p vu ■ 1 - A e-*(9^-'^) 
+-'P t 77^^ 2 5Re l——^m -i ) (8) 



1 



1-2AX + A2V //(A^> 

"• and Up = -y V. (9) 



where 



^ 1-A2 ) ^ P V1-A2 

The form ^ has the advantage that now v appears only in two places via Vp and v. In 
the superintegrable case one has (j) = (p so that ip = 1 and in (^) the second term of the 
curly bracket vanishes. Recall that on the boundary of the hermitian region (/> = or 
(/) = TT so that cos tp = ±A and we have cos ((/> — 93) = A and ip = ((1 -|- A)/(l — A))^/". 

From @ we see that for t; = 0, 00 the re-dependence (which is due only to the 
factors and Vp) disappears. For u = we get immediately 

A^,,o(. = 0,A) = 2y7^=| \^^[^ ; , (10) 

while for u = 00 the principal value integral can be calculated analytically, leading to 



= 0, vr 
0=f 



(11) 
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INVERSE TEMPERATURE A 

Figure 1: Phase diagram of the integrable Chiral Potts model showing the boundaries of the non- 
hermitian region (j> ~ 0, n (full curves). The dashed-dotted curves give the boundaries of the IC phase 

for the Z„-models with n = 3, 4, 6, 96. 

3 Boundary of the region where AEq=i becomes negative. 

As was already mentioned at the end of Sec.l, the C-IC phase transition arises where 
the ground state looses translational invariance. This is the case where AEq=i(P), by 
varying A at fixed as a function of P acquires a second order zero at some P ^ 0. 

In the regime where parity-invariance is broken, in general the second order zero 
will appear at P 7^ anyway. Therefore, even not knowing the precise relation between 
V and P, we can as well look in (^) or for such a second order zero as a function of v. 
In reft3 McCoy and Roan have used this method to obtain the IC phase boundary for 
Z3. For general Z„, since the analogue of the McCoy- Roan formula for gaps with Q >2 
has not been worked out, there are only qualitative arguments (experience from finite- 
size calculations) that it is the Q = 1-gap and not a higher-Q gap which determines 
the C-IC phase boundary. 

For the A£'Q=i-gap we did the calculation for many values of n for various (p and 
obtain the curves shown in Fig. 1. The low- and high-temperature IC-phase boundaries 
are related by duality: A ■s-^ together with (^9 

Looking at our numerical results for very large n {n = 10^ . . . 10^) we found that in 
the limit ^/ns'm{^p/2) — > 00 the phase boundary curve Xdf) and the corresponding 
value of V at the second order zero, Vc, assume the simple form (here we give the values 
Ac for the low-temperature branch): 

Ac = — sin2((^/2) + Vc = I - - cot {ip/2) + . . . . (12) 
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Figure 2; Energy- momentum plot of low-lying levels in the Q = 1-sector of the Zs-integrable Chiral 
Potts model for A = 0.7 and ip = 120° (which is close to the C-IC-phase boundary since A£'q=i 
is just getting negative at P ~ 0.27r). Dashed-dots curve: McCoy- Roan formula using our rapidity- 
momentum relation (16). Full and dotted curves: expansion to 4th and 3rd order in A, respectively. 

Symbols: finite-size data for up to = 11 sites. 



For (/? <C 1 the IC phase forms a wedge around A = 1 described by Ac— 1 = Unf'^^^^" + 
... where as = 0.0185(2); 04 = 0.0883(3); 05 = 0.153(1); ae = 0.204(2). Notice that 
for small n this is very narrow. 

Once knowing the result (12) it is not difficult to prove it analytically from (^) 
expanding A~^ and Vc — 1 in powers of q = ir/n. Recall that the n-dependence comes 
in (p|) only through and i/p. Two principal value integrals appear, which can 

be calculated explicitly. 



4 Determination of the Rapidity-Momentum Relation 

In view of the difficulties with high genus Riemann surfaces, probably it will take some 
time until someone will come up with an analytic derivation of the dependence P{v) 
for the gap AEq=i at general n, A and f. Therefore we find it worthwhile to do what 
is possible right away, which is to combine numerical and symmetry information for 
deducing the main behaviour of P{v): 

1) Numerical diagonalization of for up to = 3 ... 12 sites is easy and leads 
to a quite precise determination of the gap AEQ=i^r=oiP) over most part of the region 
< A < 0.9 and < 99 < 120°. A third order high-temperature expansion (which we 
performed up to n = 8) leads to almost the same values even up to A ~ 0.7, only at 



larger A there are problems for P < ^, see the Figures for Z3 given in ref. 0. 
All these results show that AEQ=i^r=oiP) smooth with just one maximum and one 
minimum in the first Brillouin zone of P. Also (^) for fixed n, A, (/5 as a function of 
V has only one single maximum and minimum in the whole range — cxd < v < +00. 
Therefore numerically it is easy to determine the mapping P{v) which brings both 
curves into coincidence. 

2) Applying this method to the superintegrable line, after some educated guesswork 
we find that the Z3 formula quoted in the text before eq.@, should be generalized to 
Z^i as I 

,p l-uj'^/'^v sinip ^ ^ 27r(n-l) , , 

e-'^= or v = with Poo = — ^ -■ 13) 

sini(P-Poo) n ^ ^ 

3) We want to summarize also other numerical results in a compact empirical for- 
mula (which for ip ^ ^ has to reduce to ([T3D). For this purpose we use the fact that at 
the boundaries of the hermitian region </> = 0, tt there are special reflection symmetries 
Q, which follow from eq. (B.18) of Albertini et aM, and which determine P{v) at w = ±1: 

• At the lower boundary (j) = (see Fig.l) H^'^^ is parity invariant. In (P) this is 
feature shows up as symmetry ^ under v which appears just at cf) = 0. So at 
(j) = 0, comparing also the numerical data, for < A < 1 we get 

P{v) = 2tt- P{v~^)] so that P{v = 1) = vr; P{v = -1) = 0. (14) 

• For (f) = TT the spectra have a charge-sector-dependent shifted parity symmetry 01 
which again corresponds to a symmetry of (^) against v <-> v~^. Using (B.18) of ref. 
we find at (p = n for < A < 1 

P{v) = 2tt- P{v--^) - 47r/n 
so that P{v = 1) = {n-2)ir/n; P{v = -1) = 2{n - l)Tr/n. (15) 

By making the following ansatz for Up {up is proportional to see (^), we fulfil all 
constraints (0), (^) and for = (j) = ^ obtain (|T^): 

t^P = - ''"ff with Po = -2(</. - 99)/n; P.^ = 27r - 2(0 + <^)/n. (16) 

smi(P - Poo) 



We have calculated various dispersion curves using (|g) with (|16D and compared 
these to finite-size- and perturbative data for the Z3 to Zy-models. Because of space 
limitations, we show only Fig. 2 which gives a non-trivial example. For < A < 0.4 
the agreement is so good that we conjecture (^) to be exact in the limit A ^ 0. A 
disagreement arises if we move towards the parafermionic corner A ^ 1, (/) — > 0. Here 
the momentum interval corresponding to positive v expands such that finally there is 
only a small interval around P = 27r left for the whole negative region — 00 < w < 0. 
Put differently, both Poo and Pq seem to converge to P = as we move towards A = 1 



on the curve </> = 0. The ansatz (16) does not yet take care of this behaviour. For the 
Z3-case and 0.95 < A < 0.999, (/> = a good fit can be obtained neglecting the negative 
f-region altogether using v =| tan p/^. There A.EQ=i^r=Q{P) becomes A-independent 
with a peak value of A£'q=i = 6atP = 7rorz; = l. 

''Observe that the symmetry is in v and not in v 



5 Conclusions 



Comparing the McCoy-Roan formula for the lowest gap of the integrable Z„-Chiral 
Potts model to finite-size and perturbative data and using the special symmetry proper- 
ties at the boundary of the hermitian regime, we obtain information about the hitherto 
unknown relation between the rapidity parameter and the physical momentum ojff the 
super integrable line. The numerical agreement suggests that the conjectured formula 
( p!6| ) may be exact in the high-temperature limit. We also find an analytic formula for 
the IC-phase boundary curve of the integrable Z„-Chiral Potts model valid for large 
^/nsm {(p/2). For n — > oo the whole superintegrable line will be in the IC-phase. 
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